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Abstract 

Collisional damping widths of giant dipole excitations are calculated in 
Thomas-Fermi approximation by employing the microscopic in-medium cross- 
sections of Li and Machleidt and the phenomenological Gogny force. The 
results obtained in both calculations compare well, but account for about 

25 — 35% of the observed widths in ^"^^Sn and ^^^Ph at finite temperatures. 
21.30.Fe, 24.30.Cz 
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During last several years, much work has been done to understand the damping proper- 
ties of giant resonance excitations at zero and finite temperature [0,0. In medium- weight 
and heavy nuclei at relatively low temperatures the overwhelming contribution to damping 
arises from the spreading width due to mixing of collective excitations with more compli- 
cated doorway states. There are essentially two different approaches for calculation of the 
spreading widths: (i) Coherent mechanism due to coupling with low-lying surface modes 
which provides an important mechanism for damping of giant resonance in particular at low 
temperatures f^-^, (ii) Damping due to the coupling with incoherent 2p-2h states which 
is usually referred to as the collisional damping [[7|-[T0|], and the Landau damping modified 
by two-body collisions |TlHT3|. Calculations carried out on the basis of these approaches 
are partially successful in explaining the broadening of the giant dipole resonance (GDR) 
with increasing excitation energy. In particular, the model based on the coupling with ther- 
mal shape fiuctuations in adiabatic approximation provides a reasonable description for the 
dipole strength distribution in ^osp^ ^-^^^^ ^j-^g ^"^^Sn nuclei [p!^-[l6|. In this work, we do not 



consider the coherent contribution to damping, but investigate the collisional damping at 
finite temperature due to decay of the collective state into incoherent 2p-2h excitations in 
the basis of a non-Markovian transport approach. In order to assess how much of the total 
width of giant resonance excitations is exhausted by decay into the incoherent 2p-2h states, 
we need realistic in-medium cross-sections which interpolate correctly between the free space 
and the medium. In previous investigations, incoherent contributions have been estimated 
by employing either free nucleon-nucleon cross sections or an effective Skyrme force, which 



provide, at most, a semi- quantitative description of the collisional damping widths [0. In 
this respect, the Skyrme force provides a poor approximation in the collision term, since in 
the vicinity of nuclear surface, it does not match at all to the free nucleon-nucleon cross- 
sections. In this letter, we present calculations performed in Thomas-Fermi approximation 



by employing the microscopic in-medium cross-sections calculated by Li and Machleidt |]T8 
which interpolate correctly between the free space and the medium, and provide the best 
available input for determining the magnitude of the collisional damping at finite tempera- 



ture. For comparison, we also carry out calculations using the phenomenological finite range 



Gogny force and the zero range Skyrme force |T9|,pO 



The small amplitude limit of the extended TDHF provides an appropriate basis for 
investigating collective vibrations, in which damping due the incoherent 2p-2h decay is 
included in the form of a non-Markovian collision term PTI. Then, a description for small 



density fluctuations, Sp(t) = pit) — po, is obtained by linearizing the extended TDHF theory 
around a finite temperature equilibrium state po [HiH, 

d 

^^-qI^P - [^0, 5p] - Po] = (1) 
~ j'dt'Tr2 [t;,p?p°e-'^o(*-*')[Q(t'),^^]e+''^»(*-*')(l - pl){l - pi)] - h.c. 

where 6U represents small deviations in the effective mean-field potential. The right hand 
side describes a non-Markovian quantal collision term in which the quantity Q{t) denotes 
the distortion matrix associated with the single particle density matrix, Sp(t) = [Q(t),Po]- 
It is convenient to analyze collective vibrations by expanding the small deviation Sp{t) in 
terms of finite temperature RPA modes 0\ and Ox, 

5p{t) = Y: ^A(t)[Ol,po] - zl{t)[0,,po] (2) 

A>0 

where zxit) and zl{t) denote the amplitudes associated with the RPA modes, 0\ and Ox, 
which are determined by the finite temperature RPA equation p3| , p^ . Substituting this 
expansion into eq.(l) and projecting by Ox, and noting that Q(t) = J2zx{t)0\ — zl{t)Ox, 
we find that the amplitudes of the RPA modes execute damped harmonic motion [^[T^ 

ihj^zx-hujxzx = f dt'T.x{t-t')zx{t') (3) 

where the right hand side describes the self-energy of the collective mode due to coupling 
to 2p-2h excitations. In the Hartree-Fock representation, the Fourier transform of the self- 
energy is given by. 
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Ea(u;) = — i X — r~ [nkniniUj - niUjUkrii] (4) 



4 hco — Ae + ir] 



where = 1 — n^, Ae = + — — , and i] is a. small positive number. The self-energy 
agrees with the expression written down by Landau [^], in which the real and imaginary 
parts, T,\{uj) = Aa(u;) — ^T\{uj), determine the energy shift and the damping width of the 
collective excitation, respectively |^. 

We consider the expression for the self-energy in the Thomas- Fermi approximation, which 
corresponds to a semi-classical transport description of the collective vibrations. It can be 
shown that the 2p-2h self-energy in Thomas-Fermi approximation can be deduced from the 
quantal expression by replacing the occupation numbers with the equilibrium phase-space 
density given by the Fermi-Dirac function rii — > f{e,T) = l/[exp{e — fi) /T + 1] with fi as the 
chemical potential, and summations over the 2p-2h states with integrals over phase-space, 
S — > / d^rcPpi(Pp2(Pp3(fp4: |]8|j26[|. Furthermore, spin-isospin effects can be included into the 



treatment by considering proton and neutron degrees of freedom separately. The small devi- 
ations of the density matrices Spp(t), Spn(t) of protons and neutrons are determined by two 
coupled equations analogous to eq.(l). The collision terms in these equations involve binary 
collisions between proton-proton, neutron-neutron and proton-neutron, and a summation 
over the spins of the colliding particles. Since in isoscalar/isovector modes, proton and neu- 
tron densities vibrate in-phase/out-of phase, Spp(t) = ^6pn(t), we can deduce equation of 
motions for describing isoscalar/isovector vibrations by adding and subtracting the corre- 
sponding equations for protons and neutrons. As a result, in Thomas- Fermi approximation, 
we obtain for the coUisional widths of isovector modes [Q, 

n = ^ / dp,dp2dp3dpAmpp + Wnn) (^) + 2Wpn (^^^ ]Z fM^A (5) 

and a similar expression for the isoscalar modes. Here, N\ = f drdp{xxy{ — -§^f) is a nor- 
malization, Axa = Xa(1) + Xa(2) - Xa(3) - Xa(4), Axx = Xa(1) - Xa(2) - Xa(3) + Xa(4), 
Z = [6{hLJx — Ae) — 6{huJx + Ae)]/hLJx, and Xa(^) denotes the distortion factor of the 
phase-space density Sf(t) = xx(t){—df/de) in the corresponding mode. In this expression, 
transition rates Wpp, Wnn, Wpn associated with proton-proton, neutron-neutron and proton- 
neutron collisions are given in terms of the corresponding scattering cross-sections as. 
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We apply the formula (5) to calculate the coUisional widths of the GDR excitations by 
taking distortion factors of the momentum distribution according to the scaling picture as 
Xd{p) = P ■ e, where e denotes the unit vector along the relative displacement of the proton 
and neutron momentum distributions. Due to momentum conservation, terms involving Wpp 
and Wnn drop out, and the damping is determined by the proton-neutron collision term. The 
spin averaged proton-neutron cross-section associated with an effective residual interaction 
can be expressed as. 
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where q = (pi — P2)/2, q' = (pa — P4) /2 are the relative momenta before and after a binary 
collision, and < q; S*, T|^;|q'; S*, T >a represents the fully anti-symmetric matrix element 
of the residual interaction between two particle states with total spin and isospin S and 
T. By noting that, S=T=1 and S=T=0 matrix elements of the interaction are space anti- 
symmetric, and S=l, T=0 and S=0, T=l matrix elements are space symmetric, we find 
that proton-neutron cross-section associated with the Gogny force is given. 
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where m*Q denotes the effective mass corresponding to the Gogny force, and the quantities 
If and I~ are the symmetric and anti-symmetric matrix elements of the Gaussian factor in 
the force. 



(9) 



In these expressions, p is the local density and Wi, Bi, Hi, Mi, pi denote the standart pa- 
rameters of the Gogny force |T^J2^. In a similar manner, we can calculate the cross-section 
associated with the Skyrme force and obtain [S], 
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^^^^ " ""^^^ ^ f[ai-Xo) + 1,(1 -X0(q^ + q'^) + 1(1 -X3)p"f+ (10) 



3[to(l + xo) + ^(1 + Xi)(q2 + q'') + |(1 + X3)p°]'+ 



[(|f(l-X2f + 3(^)^(1 + X,n(q-q')^) 



The expressions for the effective masses m^(r) and mj(r) for the Gogny and Skyrme forces 
can be found in [^,^]. The angle B between q and q' in the cross-sections (9) and (13) 
defines the scattering angle in the center of mass frame, and the total cross-section is found 
by an integration over this angle, crp„ = 2tt J sinQdQ (da/dQ)^^. This integral is trivial for 
the Skyrme cross-section, but a numerical integration is needed over the scattering angle to 
find the total proton-neutron cross-section for the Gogny force. 

In figures 1 and 2, we compare the microscopic in-medium cross-sections of Li and Mach- 
leidt with the cross-sections (dotted lines) of the Gogny force (dashed lines) and the Skyrme 
force (solid lines) with the SkM* parameters. For this purpose we employ a convenient 



parameterization of the microscopic calculations as reported in |18 



a"' = 131.5 + 0.092 20.2 - 1-° + °-°°34 g.^" 

pn I lab J 1.0 + 21.55 ^ ^ 

where Eiab = (pi — p2)^/2m = 2q^/m is the kinetic energy of the projectile in the rest 
frame of the target nucleon which is also equal to twice the energy available in the center 
of mass frame. Figure 1 illustrates the total cross- sect ions as a function of the bombarding 
energy Eiab at two different nuclear matter densities p = Po = 0.18/m~^ (top panel) and 
p = po/2 (bottom panel). The cross-sections shown in the left and right parts of the 
figure are calculated with the bare nucleon mass and the corresponding effective masses, 
respectively. Around the normal matter density p ^ po, and over a narrow energy interval 
around 150 MeV, these cross-sections roughly match, however the phenomenological cross- 
sections deviate strongly from the microscopic cross-sections at lower densities and lower 
and higher energies. Discrepancies are even larger in the calculations with the bare nucleon 
mass. In figure 2, the cross-sections are plotted as a function of density at the bombarding 
energy Eiab = 100 MeV. For decreasing density, the microscopic calculations approach the 
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free proton-neutron cross-section and compare well with the experimental data, where as the 
phenomenological cross-sections strongly increase and reach large values in free space. From 
these observations, we can safely state that the microscopic calculations of Li and Macleidt 
provide a more reliable description of the in-medium cross-sections than those given by the 
Gogny and Skyrme type forces. 

In a previous work, the momentum integrals in the expression (5) for the damping width 
are evaluated in an approximate manner using a method familiar from the Fermi liquid 
theory. The method subsequently improved by incorporating the surface corrections [§. 
Here, we perform the momentum integrals exactly. For this purpose it is convenient to 
transform the variable into the relative momentum and total momentum before and after 
the colhsion, q = (pi - P2)/2, P = pi + P2, and q' = (pa - P4)/2, P' = ps + P4. Then, in 
expression (5), the integral over P' is done immediately to give 
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Io{r) = J d'PdVq'iq ■ e - q' ■ e)^ (^] 5(huD ~ e' + e)hhhU (13) 

\ / pn 

where e = q^/m* and e' = q'^/m* represent energies of two particle system in the center of 
mass frame before and after the collision. In terms of the integration variables, the products 
of the phase-space factors have the following forms, 

/1/2 = 7 ^ r (14) 

l + A^{E,e)+A{E,e)[exp{+zVE^/T) + exp{-zVE^)/T) 

and 

1 + A^{E, e') + A{E, e') [exp{+z'^/E^)/T + exp{-z' ^/E^) /T) 
where A{E,e) = exp {{E + e)/2 — ep) /T, and z = cosO and z' = cosO' are the angles 
between q and P, and q' and P, respectively. In the further evaluation of the momentum 
integrals, we neglect the angular anisotropy of the cross-sections and make the replacement 
{dcr/dQ)pn (7pn/4:TT to obtain. 



rn* poo roc r + l r + l _ _ 

lDir) = {2n)\m*V^f— / ^dE^^dea^n^t + 1') / dzdz' hf2hh (16) 

6 Jo Jo J-i J-i 

where e' = e + and E = P"^ /Am* denotes the kinetic energy of the center of mass. 
In obtaining this resuh, we make the replacement (q ■ e — q' • e)^ (q^ + q'^)/3 in the 
integrand of (13), which follow from the fact that the integral of the cross terms vanishes 
and the integral is independent of the orientation of the unit vector e. It is possible to carry 
out the angular integrals over z and z' analytically, and we perform the remaining three fold 
integrations over E, e and the radial coordinate r numerically. 

In the numerical evaluations, we determine the nuclear density p(r) in Thomas-Fermi ap- 
proximation using a Wood-Saxon potential with a depth Vq = —44 MeV, thickness a = 0.67 
fm and sharp radius Rq = 1.27 A^^^ fm, calculate the position dependent chemical poten- 
tial fi{r,T) in the Fermi-Dirac function f{e,T) at each temperature, and use the formula 
hu = 80 A- for the GDR energies. In Fermi gas picture, the effective excitation 
energy of the system and temperature are related according to, 

E* = J J^sd'^d'P ' T) - /(e, T = 0)] (17) 

For temperatures small compared to the Fermi energy ep = fi{T = 0) and in sharp radius 
approximation, the relation between the excitation energy and temperature becomes, E* = 
apT'^i where ap = ATr'^/Aep is the Fermi gas level density parameter. The experimental 
temperature T* is determined from the excitation energy using a similar relation T* = 
E* /ap, where ap denotes the energy dependent empirical level density parameter 0. 
Hence, the temperature parameter in f{e,T) is related to the experimental temperature as 
T = T* ^ ap/ ap. Figure 3 shows the collisional damping width of GDR in ^"^^Sn and "^^^Ph 
as a function of the experimental temperature and comparison with data. Calculations 
performed with the cross-sections of Li and Machleidt are indicated with dotted lines. For 
comparison, we also indicate the results with the SkM* (solid lines) and the Gogny (dashed 
lines) cross-sections with the bare nucleon mass (lower panel) and the effective nucleon 
mass (upper panel), and the points are the experimental data 0. The calculations with 
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cross-sections of Li and Machleidt exhibit a weaker temperature dependence than data and 
account for about 25 — 35% of the experimental damping widths in tin and lead nuclei at 
finite temperatures. As discussed above, the calculations with Gogny and Skyrme forces are 
not reliable, since the phenomenological cross-sections do not interpolate correctly between 
the free space and the medium. Their magnitude at low densities and in the vicinity of Fermi 
energy, where the dominant contributions to damping arises, become much larger than the 
cross-sections of Li and Machleidt. As a result, the phenomenological forces predict larger 
damping, although, the magnitude of the damping is reduced by the effective mass, as 
seen from the top panel of figure 3. However, it is surprising that, the damping width 
calculated with the Gogny force including the effective mass reduction almost agree with 
the results obtained with the cross-sections of the Li and Macleidt for both tin and lead 
nuclei. Figure 4 illustrates a comparison of the calculations with the revised data taken 



from [0. In the revised data, the effective temperatures are smaller and the damping 
widths show a stronger temperature dependence than reported in 0. Consequently, the 
collisional damping accounts for a smaller fraction of the experimental damping widths. 

In order to assess the fraction of the total width of collective excitations that is exhausted 
by decay into the incoherent 2p-2h states, namely the collisional damping, we need realistic 
in-medium nucleon-nucleon cross-sections around Fermi energy. In this respect, microscopic 
in-medium cross-sections of Li and Machleidt provides the best available input, since the 
cross-sections interpolate correctly between the free space and the medium. In contrast, 
the cross-sections based on the phenomenological Skyrme force and even the finite range 
Gogny force exhibit unrealistic behavior as a function of density and energy. In the present 
investigations, by employing cross-sections of Li and Machleidt, we carry out calculations 
of the damping widths of giant dipole excitations in Thomas-Fermi approximation at finite 
temperature, and compare the results with the GDR measurements in^^^S'ra and ^osp^ nuclei. 
We also compare the results with calculations based on the Skyrme and Gogny cross-sections. 
To the extend that we accept the validity of the in- medium cross-sections of Li and Machleidt, 
we can conclude that the collisional damping of the GDR excitations is not very strong, and 



accounts for about 1/3 of F^- in both nuclei at zero and finite temperatures. Consequently, 
there is plenty of room remains for coherent damping mechanism due to coupling of the 
dipole mode with surface fluctuations. 
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Figure Captions: 



Figure 1: The proton- neutron in- medium cross-sections as a function of bombarding energy 
Eiab at several different densities. Dotted lines are cross-sections of Li and Machleidt, 
and solid and dashed lines are cross-sections associated with the SkM* and the Gogny 
forces with the bare nucleon mass (left) and the effective nucleon mass (right), respec- 
tively. 

Figure 2: The proton- neutron in- medium cross-sections as a function of density p at Eiab = 
lOOMeV. Dotted lines are cross-sections of Li and Machleidt, and solid and dashed 
lines are cross-sections with the SkM* and the Gogny forces with the bare nucleon 
mass (left) and the effective nucleon mass (right), respectively. 

Figure 3: The collisional damping width of GDR in ^'^^Sn and ^osp^ ^ function of 
temperature. Dotted lines are calculations with the cross-sections of Li and Machleidt, 
and solid and dashed lines are results with the SkM* and the Gogny cross-sections with 
the bare nucleon mass (lower panel) and the effective nucleon mass (upper panel). The 
data is taken from 01. 



Figure 4: Same as in figure 3, but shown with the revised data taken from jig . 
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